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MLWI - Malawi Land Use Survey Analysis 


DECUS Program Library Write-up 


DECUS No. 8-326 


ABSTRACT 

This program was developed to undertake the analysis of the data from the Malawi Land 
Use Survey. It calculates the proportions and areas of land in each of 15 land-use classes, 
together with their standard errors. The proportions and areas in the broad classes of "cul¬ 
tivated, " "uncultivated, " and "uncultivatable" are also computed, with their standard 
errors. 

TAPES REQUIRED 

1. Form of program tape - The program tape is written in the PDP-8 FORTRAN-D language, 
and is in the source language. There are separate program tapes for the analysis of map 
sheets, administrative areas, and natural areas. 

2. Form of data tape - The data tape should contain the identification number of the area, 
the area, and numbers of points falling in the successive land-use classes. 

OPERATING INSTRUCTIONS 

. FORT Program in high-speed reader 

*OUT-S:MLWI 

* 

*IN-R: 

■fe ^ 

T 

*READY Data tape in low-speed reader 

T 

If the program has already been compiled onto the disk, it may be called back into core as 
follows: 

. FOSL 

*IN-S:MLWI 
* 

*OPT- 
* T 

*READY 

t 


Data tape in low-speed reader 







OUTPUT 


The program prints the identification number of the area, and the area, followed by the 
number of points in each successive land-use class. These are followed by the proportions 
and areas and the standard errors for each land-use class. Finally, the proportions, areas, 
and standard errors are printed for the groups of "cultivated, " "uncultivated, " and "un- 
cultivatable" classes. 


STORAGE AND LIMITATIONS 
Normal for FORTRAN-D. 
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MLWI 




L 

C MALAWI LAND USE SURVEY ANALYSIS BY ADMIN. DISTRICTS 
DIMENSION V(16) 

20 TYPE 100 

100 FORMAT (/,/,/, "ADMIN. DISTRICT AND TOTAL AREA",/) 

ACCEPT 101,NA,TA 

101 FORMAT (I, E) 

TYPE 102 

102 FORMAT (/, "NOS OF POINTS IN SUCCESSIVE LAND USE CLASSES",/) 

sum=0.0 

DO 1 1=1,15 
ACCEPT 1,03, V(l) 

SUM=SUM+V(I) 

1 CONTINUE 

103 FORMAT(E) 

CULT=V(1 )+V (2)+V (3)+V (4) 

U NC D=V (5)+V (6)+V (7)+V (8)+V (9)+V(l 0)+V(l 1) 

UNCB=V(12)+V(l 3)+V(l 4)+V(l 5) 

TYPE 104 

104 FORMAT (/,/, "CLASS PROPORTIONS AREAS AND STANDARD ERRORS") 

DO 11 1=1,15 

P=V(I)/SUM 
A=P*TA 
Q=1.0-P 
SEP=P*Q/SUM 
IF(SEP)10,10,9 

9 SEP=SQTF(SEP) 

10 SEA=SEP*TA 

TYPE 105,1, P, A, SEP, SEA 

105 FORMAT (/, I, E,E,E,E) 

11 CONTINUE 
P=CULT/SUM 
A=P*TA 

Q=1.0-P 
SEP=P*Q/SUM 
IF(SEP) 3,3,2 

2 SEP=SQTF(SEP 

3 SEA=SEP*TA 

TYPE 106, P, A,SEP, SEA 

106 FORMAT (/,/, "CULT", E, E, E, E) 

P=UNCD/SUM 

A=P*TA 
Q=1.0-P 

SEP=P*Q/SUM 
IF(SEP) 4,4,5 
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SEP=SQTF(SEP) 

SEA=SEP*TA 

TYPE 107, P, A, SEP, SEA 

FORMAT (/, "UNCD", E, E, E, E) 

p=uncb/sum 

A=P*TA 
Q=1 .fS-P 
SEP=P*Q/SUM 
IF (SEP)6,6,7 
SEP=SQTF(SEP) 

SEA»SEP*TA 
TYPE 108, P, A, SEP, SEA 
FORMATC/, "UNCB", E, E, E, E) 
GO TO 20 



MLWI 
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MALAWI LAND USE SURVEY ANALYSIS BY MAP SHEETS 
DIMENSION V(16) 

TYPE 100 

FORMAT(/,/,/, "MAP SHEET AND TOTAL AREA",/) 

ACCEPT 101, NA,TA 
FORMAT (I, E) 

TYPE 102 

FORMAT (/, "NOS OF POINTS IN SUCCESSIVE LAND USE CLASSES" /) 
SUM=0.0 
DO 1 1=1,15 
ACCEPT 103, V(l) 

SUM=SUM+V(I) 

CONTINUE 

FORMAT(E) 

CULT=V(1 )+V(2)+V(3)+V(4) 

U N C D=V (5)+V (6)+V (7)+V (8)+V (9)+V(10)+V (11) 

U NCB=V(12)+V(l 3)+V(l 4)+V(l 5) 

TYPE 104 

FORMAT (/,/, "CLASS PROPORTIONS AREAS AND STANDARD ERRORS") 

DO 11 1=1,15 

P=V( D/SUM 

A=P*TA 

Q=1.0-P 

SEP=P*Q/SUM 

IF(SEP)10,10,9 

SEP=SQTF(SEP) 

SEA=SEP*TA 

TYPE 105,1, P, A, SEP, SEA 
FORMAT (/, I, E, E, E,E) 

CONTINUE 

p=cult/sum 

A=P*TA 
Q=1.0-P 
SEP=P*Q/SUM 
IF(SEP) 3,3,2 

SEP=SQTF(SEP) 


3 SEA=SEP*TA 

TYPE 106, P, A, SEP, SEA 
106 FORMAT (/,/, "CULT", E, E, E, E) 

P=UNC D/SUM 
A=P*TA 
Q=1.0-P 

SEP=P*Q/SUM 
IF(SEP) 4,4,5 
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SEP=SQTF(SEP) 

SEA=SEP*TA 

TYPE 1,07, P, A, SEP, SEA 

FORMAT^, "UNCD", E, E, E, E) 

p=uncb/sum 

A=P*TA 
Q=1 -0-P 
SEP=P*Q/SUM 
IF (SEP) 6,6,7 
SEP=SQTF(SEP) 

SEA=SEP*TA 
TYPE 108, P, A, SEP, SEA 
FORMAT^, "UNCB", E, E, E, E) 
GO TO 2jS 





L 

C 

MLWI 

MALAWI LAND USE SURVEY ANALYSIS BY NATURAL AREAS 

20 

DIMENSION V(16) 

TYPE 100 

100 

FORMAT (/,/,/, "NATURAL AREA AND TOTAL AREA") 

ACCEPT ljfll,NA,TA 

101 

FORMAT (1, E) 

TYPE 102 

102 

FORMAT (/, "NOS OF POINTS IN SUCCESSIVE LAND USE CLASSES",/) 

SUM=0.0 

DO 1 1=1,15 

ACCEPT 103, V(l) 

SUM=SUM+V(I) 

1 

^ 1 03 

CONTINUE 

FORMAT(E) 

CULT=V(1)+V (2)+V (3)+V (4) 

U NC D=V (5)+V (6)+V(7)+V (8)+V (9)+V (1 0)+V (11) 

104 

UNCB=V(12)+V(l 3)+V(14)+V(l 5) 

TYPE 104 

FORMAT (/,/, "CLASS PROPORTIONS AREAS AND STANDARD ERRORS") 

DO 11 1=1,15 

P=V(I)/SUM 

A=P*TA 

Q=1 .0-P 

SEP=P*Q/SUM 

IF(SEP)ljZf, 10,9 

9 

10 

SEP=SQTF(SEP) 

SEA=SEP*TA 

105 

11 

• 

TYPE 105, l,P, A, SEP, SEA 

FORMAT (/, 1, E, E, E, E) 

CONTINUE 

P=CULT/SUM 

A=P*TA 

Q=1 .0-? 

SEP=P*Q/SUM 

IF(SEP) 3,3,2 

2 

SEP=SQTF(SEP) 

3 

SEA=SEP*TA 

106 

TYPE 1J06, P, A, SEP, SEA 

FORMAT (/,/, "CULT", E, E, E, E) 

P=UNCD/SUM 

A=P*TA 

Q=1 .0-P 

SEP=P*Q/SUM 

IF(SEP) 4,4,5 
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SEP=SQTF(SEP) 

SEA=SEP*TA 

TYPE ]07, P, A,SEP, SEA 

FORMAT^, "UNCD", E, E, E, E) 

P=UNCB/SUM 

A=P*TA 

Q=1 .0-P 

SEP=P*Q/SUM 

IF (SEP) 6,6,7 

SEP=SQTF(SEP) 

SEA=SEP*TA 
TYPE ]08, P, A,SEP, SEA 
FORMAT (/, "UNCB", E, E, E, E) 
GO TO 20 



